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ABSTRACT 


The presence of H^O and CO^ in the earth’s atmosphere has for centuries 
frustrated astronomers in their attempt to "see" obj ects'' radiating outside the 
visible window. Recently, NASA has outfitted a modified C141 transport with 
a 91,5 cm reflector telescope designed to view in the infrared range from lu 
to lOOOp. The telescope is situated in a cavity which (because of the lack 
of infrared-passing windows) is operated open port. Spoilers have been 
designed which reduce turbulence- induced excitation of the cavity. 

Since the aircraft is designed to operate at altitudes up to 15 km, the 
effect of the H^O and CO^ is reduced significantly. Furthermore, the optically 
degrading influence of the large scale atmospheric turbulence on land-based 
telescopes is thus replaced by the, as yet little understood, effect of the 
turbulent shear layer resulting from the spoiler upstream of the cavity. The 
purpose of this report is to establish a mathematical model appropriate to 
describe the effect of turbulent shear layers on imaging systems as well as to 
examine the parameters of interest relevant to potential wind-tunnel experimen- 
tation. 

Because of the relative thinness of the turbulent shear layer, it is 
argued that the zeroth order geometrical optics plane wave approximation 
adequately describes the radiation field in the vicinity of the telescope. 

With rather mild assumptions, one then finds the following expressions for the 
average and mean square of the optical transfer function: 
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where 0 is defined in terms of the density correlation function, R^ , by 
P(t,u,;v,w) - Rp (v;w) - R p (t;w) - R p (u;v) + R p (t;u) 

and p= (u^+x, z^) , q=(u 1 +x,z 2 > ,r=(u 1 ,z 1 ) , s^u^z^, t=(x 1 +x,z 1 ), u=(x 1 +x J z 2 >, 
v=(x^,z^), and w=(x^,z 2 ). lu the special case that R p is Gaussian, <x> can 
be found explicitly in closed form, and a bound for <[t| 2 > obtained. 


I . INTRODUCTION 

The adverse effect of the earth’s atmosphere on man’s study of celestial 
bodies has been of concern since the first telescopes were built in the 
seventeenth century. Electromagnetic waves traverse the universe relatively 
unhindered in their travels only to be absorbed or dispersed in the last few 
kilometers of their journey to m^n’s eye. For example, the presence of H^O 
and CO^ in the atmosphere accounts for the absorption of essentially all 
infrared radiation (lu to lOOOy) except for several narrow windows in the 
spectrum up to 20 microns, 

An observatory at an altitude of 15 km would be above 95% of the H^O and 
75% of the C0 2 in the atmosphere and objects emitting in the infrared would 
be detectable. With this in mind, NASA has installed a 91.5 cm reflector 
telescope in a C-141 aircraft (ref. 1). Though not matching the height advan- 
tage of satellites, the aircraft has some advantages in cost and flexibility. 
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Since no known material is transparent to radiation at all wavelengths in the 
infrared range, the airborne telescope is designed to operate open port. This 
circumstance has created the potential for acoustic resonance of the telescope 
cavity. To avert the concomitant telescope vibrations, spoilers were placed 
just upstream of the cavity (ref. 2), and, in this configuration, the telescope 
is able to "see" in the infrared portion of the spectrum. However, the quality 
of seeing varies randomly in time due to the refractive index fluctuations 
introduced by the spoilers. Thus, the large scale turbulent fluctuations of 
the atmospheric boundary layer, which are responsible for such phenomena as 
the twinkling of stars, have been replaced by the relatively small scale 
fluctuations in the spoiler-induced free shear layer. 

The phenomenon just described is similar to optical degradation experienced 
by photo reconnaissance or earth resources missions as well as in the use of 
laser beams for communication systems. The occurrence of optical degradation 
due to "thin" turbulent boundary layers and shear layers common to so many 
varied endeavors has stimulated much work, both experimental and theoretical. 

It is the purpose of this paper to summarize and develop the requisite mathe- 
matical models which will provide a basis for an experimental investigation 
of the fundamental nature of the influence of turbulent shear layers on 
electromagnetic wave propagation, especially in the infrared range. 

What follows is a short history of relevant work performed during the 
last 20 years. The most thorough experimental investigation of the effect of 
turbulent boundary layers on light propagation has been that of Stine and 
Winovich (ref. 3) in 1956. Their results indicate that the loss in resolution 
due to turbulent boundary layers can be several times that experienced by 
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land-based telescopes resulting from atmospheric turbulence. Hufnagel and 
Stanley (ref. 4) showed theoretically that the optical degradation due to 
turbulence can be related to the distribution of the mutual intensity function 
across the aperture of the telescope. Fried (ref. 5) further distinguishes 
between M long~exposure" (corresponding to Hufnagel and Stanley's results) and 
" short-exposure" average optical transfer functions. In the present paper, 
the author will describe those aspects of the theoretical models devised thus 
far which are particularly relevant to propagation of infrared radiation 
through turbulent shear layers. 

II. THE ELECTROMAGNETIC FIELD 

The electromagnetic field is a vector field satisfying Maxwell's equations 
together with appropriate constitutive relations. (See Born and Wolf (ref. 6, 
pp, 1-3)). By assuming that the magnetic permeability is constant and that 
the dielectric constant does not vary appreciably over times of order X/c 
or distances of order X, the electric and magnetic vectors are found to satisfy 
identical wave equations. (X represents the wavelength and c the free space 
propagation speed.) Hence, each component of the field vectors satisfies the 
scalar wave equation: 

V 2 V - (n/c) 2 V = 0, (2.1) 

where V* V(x,t) represents any of the components of the electromagnetic field 
vectors, x is the position vector in the field, t represents time, and 
n« n(x, t) is the refractive index. If there is no preferential polarization 
direction, then the study of the scalar equation (2.1) suffices to yield the 
measurable quantities of interest, such as intensity or mutual coherence. 

(See Born and Wolf (ref. 6, pp. 387-392),) 
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Assuming V is polychromatic, and assuming further that V is square 


integrable, we have the following Fourier representation: 
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In the form (2.2), V is in general complex, which proves to be convenient. 
However, in the physical problem we will, of course, be interested in the 
real part of V. 

Substituting (2.2) into (2.1) we have 


n “i^t, d; 

u(x,m)e dw -(— ) 
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Assuming the required uniform convergence of the improper integrals involved, 
the order of integration and differentiation in (2*3) may be interchanged to 
give 

CO 

[V 2 u+(n/c) 2 w 2 u]e ^ Wt dw = 0, (2.4a) 

« 

w 00 

If we define k = w/c to be the free space wave number (rad/m), then 

(V 2 u+k 2 n 2 u)e * Wt dw “0. (2.4b) 

—co 

Finally, we have 

V 2 u+k 2 n 2 u = 0; (2.5) 

that is, each Fourier component of the complex signal, V, satisfies the 
Helmholtz time- independent wave equation. If the Fourier components u(x,w) 
are known, then the signal V is given by (2.2). Furthermore, if polarization 
is neglected, then intensity, I(x), at a point satisfies the relation 
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We will henceforth concentrate our attention on u(x,m) rather than V(x, t) 
since one determines the other. 

We wish now to specialize the problem to that of a plane wave propagating 
through a turbulent shear layer of thickness, L. The turbulence will be 
considered to occupy the infinite layer bounded by the planes z=0 and z=L. 

The undistorted plane wave travels in the positive z direction, entering the 
shear layer at the plane z=0 and exiting at z-L. It will furthermore be 
assumed that the wave radiates to » in the positive z direction, unaffected 
except for the turbulent region just described; in particular, it is assumed 
that the body which borders the shear layer does not reflect any of the 
incident radiation. (This is sometimes referred to as the Sommerfeld 
radiation condition.) With these boundary conditions understood, each Fourier 
component of the radiation satisfies the Helmholtz equation (2.5), 

It should be pointed out that since the refractive index n varies 
randomly in space, equation (2.5) is a random partial differential equation 
with variable coefficient, and hence it would be fruitless to seek an exact 
solution. However, in most cases of practical interest, one can write 

n(x) = l+n 1 (x), (2.6) 

where |n^|«l; i.e,, the refractive index is perturbed only slightly from 
its free space value. Because of this fortuitous circumstance, perturbation 
solutions of (2,5) abound. Most of these methods are described in references 
7 and 8. Because of the thinness of the turbulent layer under consideration. 
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however, it is the author T s opinion that diffraction effects can be ignored, 
and a geometrical optics description suffices# Evidence to support this 
contention follows in section III. 

Tatarski [8, p. 174] seeks a solution of (2.5) of the form: 

r , 1 ,1 t ike(x) 

u(x) = [u Q (x) 4* — u 1 (x) 4 ^ 2 u 2 (x) + . .#]e 

(The dependence of u on u) will henceforth be implied rather than explicitly 
stated.) He finds the following expressions for u^, and 0s 

u Q (x,y,z) - A q exp{-% 
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where V 2 denotes the lateral (perpendicular to the direction of propagation) 
Laplacian operator 
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and Aq denotes the initial undistorted wave amplitude. Tatarski argues that 
terms higher than can be neglected if (i) l/k<<l^ and (ii)/X L<<1^ , where 
1 denotes the Kolmogoroff inner scale of turbulence. Condition (i) places 
an upper limit on the wavelength 

X«2tt 1 , 
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whereas (ii) essentially restricts the thickness of the shear layer as 
follows : 


L<<l?/A. 

i 


(Restated, condition (ii) requires that the radius of the first Fresnel 
zone, /XL, be much less than the inner scale of turbulence.) In the analysis 
which follows, then, we will take: 
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(2.7) 


III. THE OPTICAL TRANSFER FUNCTION 


The question of how to describe the effect of turbulence on the quality 
of an imaging system is, to some extent, a matter of taste. Several para- 
meters have been used by designers to describe optical quality depending on 
the design aspect of interest. A description of several of these quality 
factors is given by O’Neill [9, Ch. 7]. With few exceptions, all can be 
related to the so-called optical transfer function (OTF), t, which is defined 
to be the normalized two-dimensional spatial Fourier transform of the image 
plane intensity distribution due to a monochromatic point source; viz: 
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where f and g are spatial frequencies with units, radians per meter. In 
(3.1), the function s(x,y), which denotes the intensity at the point (x,y) 
in the image plane due to a point source at (0,0) in the object plane, is 
sometimes referred to as the point spread function. If. one views an extended 
object emitting incoherent light, then 


I(f,g) = r(f,g) 0(f , g) , 


(3.2) 


where I and 0 are the two-dimensional spatial Fourier transforms of the 
intensity distributions in the image and object planes respectively. If 
one were viewing an object with randomly varying intensity distribution, 
then instead of (3.2), one has: 

<i>ii( f , g ) “ |t(f,g)| 2 <£ 00 (f,g) , (3.3) 

where <J> and which denote the Fourier transforms of the intensity 

ii DO’ 

autocorrelation functions in the image and object planes respectively, are 
the image and object Wiener spectra respectively. ($(f,g) corresponds to 
the power spectral density for time varying random functions.) To summarize, 
then, for incoherently illuminated objects, the optical transfer function 
relates the Fourier transforms of image and object intensities in a linear 
manner. 

/ 

Let us define the pupil (or aperture) function, G(x,y), to be the 
complex disturbance at the point (x,y) in the aperture plane due to a point 
source; G(x,y) will be defined to be zero if (x,y) lies outside the aperture. 
(Note that we restrict our source to lie within an isoplanatic patch of the 
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working field.) It can be shown (see, for example. Born and Wolf [6, p.485] 
or O'Neill [9, p. 77]) that the optical transfer function is related to G as 
follows : 


,kx kv. 

T( 1P R = 


G (x'+x,y '+y)G* (x 1 ,y ' ) dx'dy' 
f? G(x' ,y f )G*(x'y' ) dx'dy' 


(3.4) 


where R denotes the focal length of the system. (Note that the integrations 
in (3.4) are only formally infinite since G vanishes identically outside the 
aperture.) The representation (3.4) will be particularly convenient as we 
now proceed to describe the effect of the turbulent layer on optical quality. 

If we assume that a point source is located at a distance far removed 
from the optical system under consideration, then the radiation entering the 
turbulent layer is, to good approximation, a plane wave. But then we have: 
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(3.5) 


where 


G 0 (x,y) = 
and u(x,y,z) is given by (2.7). 

Before proceeding further, a word is in order here regarding the efficacy 
of t as a measure of optical degradation due to turbulence. By definition, 
t is the transfer function relating the intensity distribution in an extended 
object to that of the image. Its value as a transfer function is thus limited 


Jl, (x,y) in the aperture 
\ 0, (x,y) not in the aperture 
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to incoherently illuminated objects. On the other hand, if one observes a 
star with a telescope, then the intensity distribution in the focal plane 
is exactly the point spread function, whose normalized Fourier transform 
is given by (3.1). Hence, even though the use of x as a transfer function 
must be limited to the ideal, and rarely realized, case of incoherent illu- 
mination, its usefulness as a measure of optical degradation is broader. 

(The derivation of transfer functions appropriate to coherent or partially 
coherent extended objects is given in Born and Wolf [6, Sections 9.5.1 and 
10.5.3 respectively].) 

Since, by (2.7), u depends on n 1 , which varies randomly with position 
in the turbulent field, then (3.5) implies that for each fixed spatial 
frequency, x is a complex random variable. Hence, a complete description of 
x requires determination of its probability density, or, equivalently, its 
moments. In practice, one is fortunate to determine the first two moments, 
namely the mean, or expectation, and the mean square. If we denote mean by 
< >, then from (3.5), 



where 

GJ 2 ' (x'jy'jx.y) - G o (x , +x,y , +y)G 0 (x' ,y*) 

and 

M 2 (x’ ,y* ;x,y) = <u(x'+x,y f +y,L) u*(x* ,y’ ,L)> . 

An exact equality is, unfortunately, not realized since, in general, the 
mean of a quotient is not equal to the quotient of the means. It is reason- 
able, however, to accept (3.6) as a first approximation. (See [10* p. 151] 
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for the derivation of a series expansion for <t> of which (3,6) gives the 
first term.) Similarly, the mean square of t is given by: 

00 

r ff r 

kx kv JJJJ (x r ,y r ;u f ,v f ;x,y)M i+ (x t ,y* ;u* ,v* ;x,y)dx ? dy 'duMv* 

< I t( r“ , r ) l 2> * ^ (3,7) 

f G^ 4 ) (x f ,y 1 ;u’ ,v‘ ;0,0)M (x T ,y ’ ;u f ,v r ;0,0)dx , dy , du , dv' 

JJJJ 0 4 

where 

G^4) (x* ,y’ ;u’ f v' ;x,y)=G 0 (x , +x,y , +y)G 0 (x f y , )G 0 (u , +x f v , +y)G 0 (u T ,v T ) 

M 4 (x’ ,y' ;u’ ,v t ;x,y)=<u(x t +x,y , +y,L)u*(x T ,y ' ,L)u(u T ,v' ,L)u* (u'+x,v’+y,L)> . 

Note that the influence of the turbulence on <t> is felt through the second 
moments of the disturbance in the aperture plane, namely the mutual intensity 
(in the numerator) and intensity (in the denominator), whereas the stochastic 
contribution to <jxj 2 > arises from the fourth order moments of the disturbance 
in the aperture plane. Although the expectations or means which occur here 
are ensemble averages, they may be calculated in practice by taking time 
averages if the light disturbance can be assumed stationary, an assumption 
almost surely warranted. 

Equations (3.6) and (3.7) are quite general, and do not depend on the 
particular model of the disturbance. If we utilize the geometrical optics 
approximation, described in section II, then the mutual intensity is given by: 

M 2 (x , ,y';x,y) = | A Q | 2 <e A+lB >, (3.8) 
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where 


A - XCx'+Xjy^jL) + X(x* ,y ’ ,L) 
fL 


- k ! (L-?) [v|n 1 (x'+x^y’+y ,0 + V^n 1 (x’ ,y ' , C) ]dC 
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B - k[Y (x , +x,y , +y > L) - Y(x\y r ,L)] 

,L 


[n^x'+xjy'-ty, c) - n x (x* ,y ' , 0 ] d£. 


Similarly, the fourth order correlation function occurring in (3.7) is 
given by: 


M 4 (x' .y'ju'.v'jx.y) = |aJ 4 <e A +lB >, 


(3.9) 


where 


A’- -h 


(L-C) [V^n^ (x f +x,y f +y ,C)+V|n 1 (x 1 ,y' , ^;)+v|n 1 (u , ,v f » C)+v£n x (u'+x,v T +y, 0 ]dC 


and 


B r =k 


[n 1 (x^Xjy’+y, C)-n 1 (x* ,y T ,C)+n 1 (u’ ,v' , 0-n 2 (u f +x,v f +y, 0 ]d5. 


Unless we know the probability density for n^, we have reached an impasse. 
However, if we are willing to assume that the variables A, B, A 1 , and B* are 
Gaussian, then further progress can be made. In particular, we will assume 
that A and B are jointly Gaussian and similarly for A* and B’. (From the 
definitions of A, B, A T , and B* , an appeal to the central limit theorem 
seems to lend credence to such an assumption.) Then, by carrying out the 
integration inferred by < > with respect to the Gaussian density, one finds 


<e A+iB > B e <A>+i<B> e %(c^ - a h e ^ cov(A,B) 


(3,10a) 
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and, similarly. 


A'+iB* <A’>4i<B , > %(a^ t - a gi) i covCA^B*) 
<e > s= e e e 


(3.10b) 


where o 2 denotes variance and cov denotes covariance. In order to carry out 
the calculations indicated by (3.10), we proceed to determine the joint 
moments for A and B, (Exactly similar results will follow for A* and B f .) 
For convenience in notation, we make the definitions: 


»y ^ » z ^ jx^ ,y ^ » z^) = 

<n 1 (x 1 ,y 1 ,z 1 )n 1 (x 2 ,y 2 ,z 2 )> 

(3.11a) 

P(x 1 ,y l ,z 1 ;x 2 ,y 2 ,z 2 ) = 

<vfn x (x^y^z^vlnj (x 2 ,y 2 ,z 2 )> 

(3.11b) 

Q(xi,y 1 ,z 1 ;x 2 ,y 2 »Z2^ ” 

< [vjnj (x^,y^,Zj)]n^ (x 2 ,y 2 » z 2 )> • 

(3.11c) 


Because the operation of ensemble averaging is commutative with both inte- 
gration and differentiation, then from (3.8) 


<A> = 


(L-?) [V|<n 1 (x'+x,y T +y, 5)>+v|<n x (x f ,y 1 ,C)>3^C 
<B> ■ k | [<n 1 (x , +x,y'+y»C) > - < n 1 (x' ,y T , ?)>]dC . 

rw 
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fL 


( 3 . 12 ) 


Note that if one assumes <n^> = constant throughout the turbulent field 
(a not unreasonable assumption), then <A> ** <B> = 0. Furthermore, from 
(3.8) and (3.11) one finds: 
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[R(x f .y'.CjSx'.y' ,C 2 )-R(x'+x > y 1 +y, C 1 ;x , ,y , ,? 2 ) 
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(3.13b) 


and 

<AB> =» k/2 


LL 

’ 
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CL-t;^) [Q(x* ,y 1 , ;x* ,y* , c; 2 ) 

+Q(x’+x,y ? +y, ;x* ,y* , 5 2 )-Q(x’ ,y T , ;x'+x,y f +y , C 2 ) 

«Q(x'+x,y ? 4y, ? 1 ;x l +x,y , +y,C 2 )3 d^d^. 


(3.13c) 


Similar expressions hold for A f and B*. Clearly, without further simplifying 
assumptions the problem remains, in a practical sense, intractable. To 
elaborate, observe that from (3.11) 


P(x 1 ,y 1 ,z 1 ;x 2 , y 2 , z 2 )= 


Q(x 1 ,y 1 ,z 1 ;x 2 ,y 2 ,z 2 ) - 


3 4 R 3 4 R 3 4 R + 3 4 R 

3x 2 3x| 3yfsx 2 3x 2 3y 2 3yf9y| 

3 2 R _9^R 

3x^ 3y2 


(3.14) 


Hence, even if one were able to experimentally determine the refractive 
index (density) correlation function, evaluation of <A 2 > and <AB> would 
require a numerical differentiation of fourth order and second order respec- 
tively. Although this is feasible in principle, it is very difficult to do 
with any accuracy in practice, without a very fine mesh measurement scheme. 

To conclude our general discussion, it should be noted that in practice 
(see, for example, [11]), it is the modulus of the OTF, called the modulation 
transfer function (MTF), which is measured, rather than the OTF itself. The 
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statistics of the MTF are related to those of the OTF as follows: 

0 < | <t> ] < <]x|> < 1. (3.15) 

Hence, if the expressions for < t> obtained in this section are used to 
estimate the average MTF, the estimate can be expected to be con- 

servative. Furthermore, the variance of the MTF is bounded as follows: 

0 < a 2 | = <|t| 2 > - <|t |> 2 £ < | t | 2 > -• |<t >[ 2 «* a 2 . (3.16) 

Knowledge of the variance, together with application of Chebyshev’ s 
Inequality (see, for example, [12, p. 20]), yields information about the 
variability of "seeing” conditions. In particular, 

Pr(|t - <t> |> aa^) < 1/a 2 

and 

Pr( ] t | - <]t|>|> a<? | T |^ 5 1/w 2 » 

for the OTF and MTF respectively. 

In concluding this discussion, it is interesting to note that from 
(3,15), the MTF suffers less degradation on the average than the OTF. 
Furthermore, from (3,16), the variance of the MTF is, in general, smaller 
than that of the OTF, 

Homogeneous Turbulence . If one is willing to assume that the refractive 
index spatial variations are statistically homogeneous, then certain simpli- 
fications can be made. In particular, <n x > * constant and, hence, from 
(3.9) and (3.12), 

<A> 5:5 <B> - <A T > = <B’> = 0. 
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Furthermore, the correlation function R is now a function only of the 
difference in the coordinates rather than the actual coordinates of the two 
points under consideration; i.e.. 


R(x 1 ,y 1 ,z 1 ;x 2 ,y 2 ,z 2 ) = R(x 2 «x 1 .y^y^z^z^ . 


(3*17) 


Then, the relations (3.13) become: 

L 

<A 2 > - k | (L-C 1 )(L-C 2 )[2F(0,0;c 2 -C 1 )+P(x,y,C 2 “? 1 ) 

0 

+P(-x,-y,C 2 -C 1 )] d? x dC 2 


(3* 18a) 


<B 2 > = k 2 


[2R(0,0,C 2 “C 1 )“R(x,y,C 2 -e 1 )“R(-x,-y i ,C 2 -C 1 ]dc 1 d? 2 (3.18b) 


<AB> = (k/2) 


(L-C^fQ^x.-y.C^C^-QCxjy,?^?^] dJ^d^, (3.18c) 


where, if R=R(£,n, 0 » then 


3 2 R , „ 3 4 R , 3 4 R 

Ptt.n.C) = + 2 + 

3£ 4 3C 2 Sn 2 


Q(£,n,0 = — + — 

3C 2 3n 2 


(3.19) 


Hufnagel and Stanley [4, Theorem II, p. 61] have shown that 
Q(-x,-y, ^ 2 -^ x ) = Q(x,y,C 2 -C 1 ), 

which implies that <AB> - 0 and, hence, that A and B are uncorrelated. Some 
simplification in the expressions for <A 2 > and <B 2 > can be achieved by making 
a change of variable (see Appendix). In particular, the double integrals 
are reduced to single integrals as follows: 
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<a 2 > 

<b 2 > 


L 

(1/12) | (2L+z)(L-z) 2 [2F(0,0 t z)+P(x,y,z)+P(x,y,-z)3dz 

0 

L 



(L-z) [2R(0, 0, z)-R(x,y ,z)-R(x,y ,-z) ] 


dz . 


(3.20) 


If we should assume further that isotropic conditions prevail, then 
R(x,y , z)=R(x,y,-z) and P (x,y , z)=P (x,y ,-z) and (3.20) can be simplified even' 
further. (In fact, one need only assume that the correlation function is 
an even function of z to obtain the same result.) 

Similarly, from (3.9), one finds: 


L 

<A f 2 > - 2<A 2 > + (1/12) | (2L+z) (L-z) 2 {2[P(u , -x T ,v , -y , ,z) 

+P(u , -x I ,v T -y ’ ,-z) ] +P(u 1 -x f -x,v 1 -y , -y,z) 

(3.21a) 

+P (u T -x f -x, v*-y *-y ,-z) +P (u’-x'+x, v T -y*fy, z ) 
+P(u t -x f +x,v , -y , +y,-z) } dz 


and 
<B’ 2 > 


2<B 2 > - 2k 2 


(L-z){2[R(u , -x f ,v $ -y’ ,z)+R(u f -x f ,v'-y' ,-z)] 


-R(u f -x , -x,v , -y , -y,z)-R(u , -x , “X,v , -y t -y,-z) 
-R(u ? -x , +x,v f -y f +y, z)-RCu t -x , +x,v T -y r -ty,-z) } dz 


(3.21b) 


Finally, combining (3.6), (3.7), (3.8), (3.9), (3.10), (3*20), and 
(3.21) we have: 


<T> = T T , 

0 R 


(3.22) 


where 


T R =exp{-% I ( 2L+z ) Q z. z l. f 1 ( x ,y,z)+2k 2 (L-z)G 1 (x,y,z)jdz} 
F 1 (x,y,z) * 2P (0,0, z) - P(x,y,z) - P(x,y,-z) 


G 1 (x,y,z) = 2R(0,0,z) - R(x,y,z) - R(x,y,-z) 
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and Tq denotes the diffraction limited transfer function of the imaging 


system. The function t , representing the random influence of the turbu- 
lence, is independent of properties of the imaging system. This decomposition 
of <t> into the product of the diffraction limited OTF and a turbulent de- 
gradation is, in general, possible only if an assumption like homogeneity 
(or local homogeneity) is invoked. Furthermore, 

CO 

|| | (x f ,y' ;u’ ,v r ;x,y)e U ^ X ,y ,U ,v ,X,y,L ^dx’ dy’du'dv 1 


<]t| 2 >= t 2 A 2 


, (3.23a) 




f G^ 4) (x , ,y l ;u f ,v r ;0,0)e U(x,,y,,uf ^'^^'^dx’dy'du’dv' 


where 


U(x* ,y 1 ,u' ,v f ;x,y,L)=% { [<A* 2 >-<B' 2 >]- 2[<A 2 >-<B 2 >]}, 


and A denotes the aperture area. If we invoke the second mean value theorem 
P 

for integrals, then (3.23a) can be rewritten as follows; 

L 

... * 

< | T |2 >BT 2 T 2 exp |^| {(2 L -<- z )( L -z) [F 2 (s 1>tl ;x t y,z)-F 2 (s 2 ,t 2 ;0,0,z) j 


+2k 2 (L-z)G 2 (s 1 , t 1 ;x,y,z) }dz V , 


(3.23b) 


where 


and 


F 2 (s,t;x,y,z) « P(s-x,t-y,z) + P (s-x, t-y,-z) + P(s+x,t+y,z) 
+P (s+x, t+y,-z) -2[P(s,t,z) +P(s,t,-z)] 


G 2 (s,t;x,y,z) = 2[R(s,t,z) + R(s,t,-z)] -R(s-x,t-y, z) 


-R(s-x,t-y,-z) -R(s+x, t+y, z) -R(s+x, t+y ,-z) . 
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Since the mean value theorem is nonconstructive, we can guarantee the 
existence of s^ and t_^ (i=l,2) which validate (3.23b), but no algorithm 
exists for determining their values. For circular and rectangular apertures, 
one can argue, however, that 


and 


0 < sj + t 2 £ (D-r) (D+r) 
i 1 



D ) 


0 < s 2 + t 2 < 

1 1 


(D -x) 2 +(D -y) 2 
x y 


< 


D 

x 



9 


respectively, where r 2ss x 2 +y 2 . Clearly, equation (3.23b) is suggestive, 
but not computationally useful. On the other hand, in the following 
example, (3.23b) will be helpful in constructing an interesting upper 
bound for < | t j 2 >. 

Homogeneous Turbulence with Gaussian Correlation , Suppose the 
refractive index correlation function is Gaussian; i#e.. 


R(x,y,z)= : 


o*- e 
n i 


] 

+< n ^> 2 . 


(A brief discussion of this example can be found in [13, p • 46 3 * ) Then 
from the definition (3.19) of P, one finds 


2 ^ 2 2 

P(x,y,z)= al {(4) [4(f) -12(f) +3] + 2 A) (£ 2 ) [1-2 (f) ][l-2(f) ] 

Tl^ a a a cl D <i D 


+ (^ 2 ) 2 [4 (f) l *-12(f) 2 + 3 ] } exp{-[(f) 2 +(f) 2 +(f) 2 ]}. 
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It might be noted that unless a=b=c, this model is homogeneous, but not 
isotropic. Even so, the correlation function is even in its arguments and 
so R(x > y,z)=R(x,y,-z) and P (x,y ,z)=P(x,y,-z) . Then, calculation of the 
constituents of as given by (3,22) gives: 


1 L ,1 


2 , ,2 2 


F 1 (x,y,z)dz --y (3(p>) +2(| 2 )(~ 2 )+3(p>) 


2^ x 4 „2 22 v 2 v 2 

+{(f 2 > [4(f) -12(f) +3] + 2(f 2 )(f 2 )[l-2(f) ] U-2^ 3 


(3.24) 


+ ( f 2 ) 2 [4(f) 4 - 1 2 ( f ) ‘+ 3] } exp{-[ (f) " + (f)"]}/I z , 


2 „ 2 


I * (2L+z) (L-z) 2 e'^'^dz 

Z 0 J 


=3c 4 {(/ir73) (L/c) 3 erf(L/c)+(l/3)[(L/c) 2 -%]e' (L/c) +(l/6)-%(L/c) 2 ). 

(NOTE; erf denotes the error function, which is defined in the usual way as: 


erf (x) = (2/ /rr) e" U du. ) 


Similarly, 


L 2 2 

2k 2 (L-z)G 1 (x,y,z)dz « 4k 2 a 2 [l-exp{-[(“) + (^) ]}] , (3.25) 

i* 1 


where 


J = f (L- Z )e" (z/C)2 dz 
2 0 J 


(c z /2){>^r(L/c)erf (L/c) - 1 + e } , 
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In the special case a=b, then, from (3.22), (3.24), and (3.25), we find 


T R =exp [-0^(4 (f ) " t 2 (e“ <r /a) 2 -l)+(r> 2 1 (~) -41 - r/a)2 ]I,/3c“ 


'a a 


+2k 2 c 2 (l-e“ (r ^ a ^ 2 ) J /c 2 } 1 . 


(3.26) 


Several calculations of <x> for this case were carried out on a CDC 3300, and 
the results are illustrated in figures 1 through 3. Some tentative conclusions 
can be drawn. In particular, high degradation can be attributed to: 


(i) large diameter optics, 

(ii) low altitudes, 

(iii) thick shear layers, 


(iv) short wavelengths. 

As might have been anticipated, the contribution of the F^ term appearing 

in the integral in (3.22) (or the I term in (3.26)) is negligible. The term 

z 

containing the wave number, k, which is very large in the geometrical optics 
approximation, essentially "swamps" the result. It is to be expected that, 
similarly, the F 2 term occurring in (3.23b) will prove negligible when com- 
pared with the k 2 term. If one ignores the contribution of the F 2 term in 
(3.23b), it is easily seen that 


<|t| 2 > < t^t 2 exp{4k 2 c 2 o 2 ^ [1-exp {-[ (—) +(£■) ]} J z /c 2 , (3.27a) 

and, hence, in the special case a=b , 

< 1 T | 2 > < t*. (3.27b) 
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Finally, for the case a=b, from (3,16), (3.22), and C3,27h), we have 


0 * ff jr| * °t < T 0 " t 0 t £ = t 0 (1 - t R> < T 0- 

It follows that 


(3,28) 


3 a. | < 3 t (3.29) 

It is the bound (3.29) which gives the 3a boundaries of figure 4. Equation 

(3,29) leads to the interesting conclusion that the smaller the influence of 

the turbulence (i » 1) on the average MTF, the smaller the variance of the 
R 

MTF. Indeed, in addition to the data presented in figures 1 to 4, the author 
performed calculations for propagation at A = lOOy for all permutations of 
the parameters D- .01 and .915 meters, a-b=2c ~ ,007 and .2 meters, L *= ,1 
and .3 meters, and <p> = .905 kg/m 3 (3km) and ,187 kg/m 3 (15km). The influence 
of the turbulence in these configurations was found to be negligible, and so 
the optical system is diffraction limited within small fluctuations. Finally, 
note that figures 1 through 4 are conservative in the sense that the correct 
MTF would, in general, be larger and the true deviations smaller than 
illustrated. 

The influence of scale of turbulence is somewhat more complicated than 
the effects described above. In particular, an asymptotic analysis of the 
expression containing k 2 on the right hand side of (3,26) leads to the follow- 
ing representations of t i 

K 
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exp(-B) ; 


L/c«l, D/a>>l, 
L/c, D/a«l, 


x R - exp{-£(r/D) 2 (D/a) 2 }*, 

x n ~ exp{-J3(L/c)“* [/rT - (L/c) - *]}; L/c, D/a>>l, 

R 

t d ' exp{-g(r/D) 2 (D/a) 2 (L/c) -1 [*^r - (L/c)" 1 ]); L/c»l, D/a«l, 
K 


where 

6 = [2tf ct u (L/X)] 2 . 

1 

If $>>1, then clearly the first case cited experiences the worst degradation. 
That is, if the upstream spoiler induces a shear layer which exhibits rela- 
tively low frequency fluctuations in the propagation direction and, at the 
same time, relatively high frequency fluctuations in the plane perpendicular 
to propagation, then degradation can be quite large. On the other hand, the 
converse situation (i.e«, L/c>>l and D/a<<l) tends to mitigate against large 
optical degradation. 

Locally Homogeneous Turbulence . It is, of course, rare that turbulence 
satisfies exactly the homogeneous hypothesis. Perhaps the simplest variety 
of nonhomo geneous turbulence is that labeled locally homogeneous. In this 
case, it is assumed that while the refractive index field is not necessarily 
homogeneous, the spatial increments of n^ are. The relevant correlation 


function is then 



which can be rewritten as follows: 


P( x 1 > x 2 ;x 3 ,x 4 ) =» %{D(x 1 ,x 4 )+DCx 2 ,x 3 )-D(x 1 > x 2 )-D(x 3 ,x 4 )} > 


(3.31) 


where 

D(^i*^j > = < [ n i(Xj) -n i(x i )] 2> 


is called the structure function of the random field, n-^. If the field 
is locally homogeneous, then 


D(x 


: ,x.) = D(x -x. ,y.-y . ,z.-z.) 
i y j i 3 i* J i 


(3.32) 


Our purpose here is to calculate the statistics for the OTF in terms of the 
structure function. 

Based on the results for homogeneous turbulence discussed above, it 
seems reasonable to assume <n]_> K constant, and, hence, that 

<A> - <B> - <A f > “ <B*> = 0. 


Furthermore, we will ignore the (presumably small) contributions of <AB>, 
<A’B'>, <A 2 >, and <A ,2 > to the statistics of t. Hence, we need only derive 
expressions for <B 2 > and <B l2 >. But, from (3.8), 


<B 2 > - k 2 


P(x f ,y 1 ,C 1 ;x T ,y ? ,^ 2 ;x*+x,y , +y,^ 1 ;x , +x,y ? +y,^ 2 )dc 1 dc 2 


(3.33) 


L, 

%k 2 Jj tD(x,y,C 2 -C 1 )+I»<-x,-y,i; 2 -S 1 )-2D(0,0,C 2 -? 1 )]dC 1 dC 2 . 


By a change of variables (see Appendix), we can write: 
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<B 2 > = k 4 


(L-z) [D(x,y,z)+D(x,y,-z)-2D(0,0,z)]dz. 


(3.34) 


Similarly, one finds: 


L> 

<B f2 > = 2<B 2 > + k 2 J (L-z)^2{D(u'-x r jv'-y 1 ,z)+D(u l -x r jv’-y 1 ,-z)] 


HDCu'-x’-XjV'-y *-y, z)-D(u , “X , -x, v’-y T -y,-z) (3.35) 


-DCu'-x^XjV'-y’+y , z)-D(u , -x , +x,v 1 -y '+y ,-z)? dz. 


Finally, combining equations (3. 6)- (3. 10) , (3.34), and (3.35), we have: 


<T> = Vr * 


(3.36) 


where 


r ■ exp{-k 2 /2 (L-z) [D(x,y,z)+D(x,y,-z)-2D(0,0,z) ]dz); 


R 


also, 


C | x I 2 > = (t 2 /A 2 ) [ J (x* ,y t ;u f ,v T ;x,y) e V ^ x * u ,V ^^^^dxMy 1 du'dv' 

R p JJJJ 0 


(3.37a) 


where 


V(x* ,y ' ,u* ,v* ;x,y,L) = -%[<B ?2 > - 2<B 2 >]. 


If we again invoke the second mean value theorem for integrals, then 


< t 2 > = t 2 t 2 exp{-k 2 /2 
OR 


(L-z) H(s,t;x,y,z) dz} , 


(3.37b) 


where 


H(s, t ;x,y,z)=*2[D(s,t,z) + D(s,t,-z)] - D(s-x,t-y,z) 

-D(s-x, t-y,-z) -D(s+x,t+y,z) -D(s+x,t47,-z) . 
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The remarks following equation (3.23b) concerning s and t are applicable 


here. 

If the turbulence is assumed to be locally isotropic as well as locally 
homogeneous, then (3.34) and (3.35) become 


<B 2 > “ 2k 2 


(L-z) [D(x,y ,z) -D(0,0, z) ] dz 


and 


<B i2 > = 2{<B 2 > + k : 


0) 


(L-z) [2D(u'-x' ,v f -y 1 , z) 


-DOi'-x'-x^'-y'-yjz) -D(u , -x t +x,v , -y f +y,z)]dz), 


respectively. Hence, the above results for '^t> and <[t| 2 > simplify to: 

L 

<t> = t t « t a exp{-k 2 (L-z) [D(x,y,z) -D(0,0,z)] dz) 

U K u 


and 


< 1 T i 2 > 


(t r /k P j 


Gp 4 hx',y';u' ,v';x,y) 


X exp{-k 2 


(L-z) [2D(u , -x I ,v'-y f ,z) -DCu’-x'-XjV'-y '-y, z) 
-D(u , -x'4-x,v 1 -y r -ty, z) 3 dx f dy , du’dv t } 


t 2 t 2 exp{-k 2 


(L-z) [2D(s,t,z) -D(s-x, t-y,z) -D(s+x, t+y, z) ]dz>. 


A few remarks are in order here concerning the relation between correla- 
tion functions and structure functions. If the correlation function is known 
then the structure function is given by: 


Unfortunately, there is no analogue for recovering the correlation function 
from a measured structure function (unless the turbulence is homogeneous, 
in which case the two statistics are equivalent). Hence, in general, it is 
recommended that the correlation function be measured, rather than the 
structure function. The structure function can be calculated using (3.38) 
if it should prove of interest. If the sole object is to verify the simpli- 
fied model presented here, then either correlation or structure functions are 
acceptable, since the structure function occurs naturally in <B 2 > and <B ,2 > 
(see equations (3.34) and (3.35)). A thorough discussion of structure 
functions and their relation to correlation functions can be found in 
[14, p. 86] and [8, pp. 13-38], 


IV. CONCLUDING REMARKS 

Several conclusions and recommendations can be drawn, and observations 
made, based on the above analysis: 

1. If one assumes from. the start that 


and 


(i) <n i > “ constant. 


(ii) <A 2 >, <A i2 >, <AB>, and <A’B 1 > « 1, 


then the plane wave solution of the geometrical optics approximation yields: 


<X (^|*^) >SS ( 1 /A p ) GQte'+^y'+yjGQte’ ,y ' ) exp{-k 2 /2 


P(t,u;v,w)d? 1 d? 2 } dx * dy * , 


(4.1) 


where t-te'+x^+y, t^) , u=(x'+x,y , +y, £ 2 ), v^Cx 1 ,y' , ^) 9 and w=(x* ,y* , ? 2 ) . 

From the definition, (3.30), V is clearly related to the correlation function 


of as follows: 
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P(t,u;v,w) « R(v;w) - R(tiw) - R(u;v)+ R(t;u). 


Similarly, the mean square < x z > is given by 


< | t ( J f, i f)| 2 > = (1/A 2 ) JJJj G-'(x’,y';u',v'; x ,y) 


J 0 


X exp{-k z /2 


[P(t,u;v,w)+ P(t,s;v,q) 


(4.2) 


+ V (r,u;p,w)+ P(r,s;p,q)]d^ 1 dC 2 dx’ dy 1 du f dv 1 , 


where t, u, v, and w are defined as above, and p= (u'+x, v’+y , C j ) > 
q=(u , +x,v , +y,C 2 ) , r= (u f ,v' , ? 1 ) , and s- (u ? ,v' , C 2 ) * Hence, by measuring 
correlation functions for n^ (or, equivalently, density) on an appropriate 
mesh, one can evaluate <x> and < | x | z > by the quadratures (4.1) and (4.2). 

If the density fluctuations are homogeneous or locally homogeneous, then, 
of course, the simpler relations (3.22)-(3.23) or (3.36)-(3, 37) , respectively, 
hold. It might be noted that the more realistic spherical wave approximations 
could be expected to imply less serious degradation than indicated by the 
plane wave theory discussed here. (See [4, figure 8, p, 60].) 

2, Since n is related to the density of the shear layer by 

n 1 « Kp , 

where K~. 000223 m^/kg is the Glads tone-Dale constant, then 

R(x.’,x ( ) ® K 2 R (x.;x.), 
i 3 P i 3 

Hence, all of the results discussed above can be written in terms of density 
correlations instead of correlations of n^. The measurement of the statistics 
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of p is a nontrivial endeavor, and much activity has been devoted toward 
that end. Hot wire anemometers [15], laser velocimeters [16], crossed beam 
techniques ([17], [18], [19]) and sphere probes [20] all offer (either 
separately or jointly) some hope for carrying out such measurements. 

3. As noted earlier (see (3.6) and (3.7)), determination of the mutual 
intensity, Mr,, and fourth order moments, M^, of the light disturbance in 

the aperture plane yields the pertinent statistics of t. Hence, the moments 
of the disturbance provide an even more fundamental description of optical 
degradation than does the OTF. Some effort has been expended towards the 
measurement of such statistics. Typical results are described in [21], [22], 
[23], [24], and [25]. It might be noted that Kelsall f s shearing interfero- 
meter [11] is capable of measuring the mutual intensity, M £ , under certain 
conditions. 

4. It is interesting to conjecture that if the influence of the turbulence 
on the optical performance were known, then the resulting image might be 
"corrected" and the turbulence influence eliminated. In fact, research has 
been carried out to just this end C[26], [27], and [28]). Korff [27] makes 
the argument that <|t| 2 > may provide more useful information than does <T> 
in that signal/noise is reduced for high spatial frequencies. 

5. As remarked earlier, the choice of the OTF, t, as the indicator of 
optical quality is somewhat arbitrary. In fact, if one were concerned with 
propagating, rather than imaging, systems, the Strehl intensity, S, is perhaps 
more appropriate. It is defined to be the ratio of the maximum intensity (in 
a particular plane of observation) to the maximum intensity which would have 
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been obtained had no disturbance been present. (See, for example. Born 
and Wolf [6, p. 462] and O’Neill [9, p. 106]). In terms of OTF, we have 


S 


jji(f,g) df dg 

_DO 



(f,g) df dg 


(The definition of S is thus the analog of system bandwidth for time-varying 
systems.) Hence, if one can calculate T and its moments, similar results 
can be obtained for S. In particular, 


< S> 


00 

/ 

< T(f,g)> df dg 

iJ 


f T (f,g) df dg 


and 


<| s| 2 > 



<t(f ,g)x*(f ' ,g’)> df dg df’ dg’ 


00 



—00 


T 0 (f,g)T 0 (f ' ,g’) df dg df 1 dg f 


Note that <|s( 2 > is not directly related to <|t| 2 > . 

6. Finally, the evaluation of experimental results and comparison with 
full scale results requires knowledge of scaling laws. A derivation of such 
similarity laws is given in [13, p. 100]. In addition to the usual aero- 
dynamic parameters such as Mach number and Reynolds number, it is clear from 
the analysis of part III that the ratio of model characteristic length to 
the diameter of the receiving optics should be maintained. 
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V. APPENDIX 


We wish to derive the following formulae: 


LL 


F(y-x) dx dy 


(L-O [F(5) + F(-03 d£ 


(A.l) 


00 


LL 

t 




(L-x) (L-y)G(y-x)dx dy « (1/6) 


<2L+C)(L-e» 2 [G(£)+G(-C)]dC. (A. 2) 


To this end, we make the change of variables 


K * y-x, n - L- (x+y ) / 2 . 


(A. 3) 


Then, in the £-n plane, the double integrals (A.l) and (A. 2) are performed 
over the area indicated in figure Al. 


K 



FIGURE Al 


32 



From (A* 3), we have 


x- L-n - C/2, y= L-n + 5/2 ; 


Since the Jacobian J = 3 (x,y)/ 3 (5,rj) ■ 1, then the double integrals of 

(A,l) and (A. 2) can be written: 

0 L+5/2 L L-5/2 

F(y-x)dx dy * | F(C) [ j dn ] d5+ j F(5) [ J dn ] d£ 

00 -L -5/2 0 5/2 


LL 

r 


(A. 4) 


L L-C/2 

[F(?)+F(-C)] [ J dn J d5 
5/2 


(L-5)[F(5)+F<-5)] d5 


and 


LI 


00 


(L-x) (L-y)G(y-x)dx dy 


6(5) 


-L 


L+5/2 

(ri+5/2)(n-C/2)dn J d5 
C/2 


- L-5/2 

G(5) [ j (n+5/2)(n-5/2)dn ] d? 
3 -5/2 


(A. 5) 


L-5/2 


| [G(5)+GC-5)3 [ J(rr+5/2)(ri-5/2)dnJ d5 . 

o 5/2 


But L- ^ 2 

(n+5/2)(n-5/2)dn 

5/2 


(1/6) (2L+5) (L-5) 2 9 and, from (A. 4) and (A. 5), 


formulae (A.l) and (A. 2) follow immediately. 
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